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Ioannis  Kanellakopoulos,  Petar  V.  Kokotovic  and  Riccardo  Marino. 


ABSTRACT 

A  new  direct  adaptive  regulation  scheme  is  proposed  for  nonlinear  systems  satisfying  extended 
matching  conditions.  Two  sets  of  additional  conditions  are  given:  first,  for  the  regulation  to  be  glo¬ 
bal,  and  second,  for  asymptotic  parameter  convergence.  It  is  also  shown  that  the  new  scheme  is 
robust  with  respect  to  unmodeled  dynamics. 
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1.  INTRODUCTION 

In  the  last  couple  of  years  an  adaptive  nonlinear  methodology  has  emerged  which  combines  basic 
concepts  of  feedback  linearization  (Jakubczyk  and  Responded  1980,  Hunt  et  al.,  1983,  Isidori, 
1985)  with  the  Lyapunov  design  of  adaptive  linear  control  (Narendra  and  Annaswamy,  1989,  Sas- 
try  and  Bodson,  1989).  This  methodology  assumes,  first,  a  linear  parametrization  (the  unknown 
constant  parameters  appear  linearly  in  the  nonlinear  model  of  the  plant  to  be  controlled),  and, 
second,  a  form  of  feedback  linearization  (complete  or  partial)  for  all  values  of  the  unknown  parame¬ 
ters.  The  results  obtained  thus  far  make  two  types  of  additional  assumptions  under  which  specific 
control  laws  are  designed.  The  first  type  of  assumptions  further  restricts  the  form  of  the  depen¬ 
dence  of  plant  models  on  unknown  parameters:  in  Taylor  et  al.  (1989)  the  model  satisfies  a  strict 
matching  condition,  while  in  Nam  and  Arapostathis  (1988)  it  appears  in  the  so-called  pure-feedback 
form.  The  second  type  of  assumptions,  made  by  Sastry  and  Isidori  (1989)  and  Sastry  and  Kokoto- 
vic  (1990),  do  not  impose  further  restrictions  on  parameter  dependence,  but,  instead,  restrict  the 
nonlinearities,  assuming  that  they  satisfy  a  “linear  growth”  condition.  The  adaptive  control  laws 
designed  under  these  assumptions  differ  in  their  applicability  and  complexity.  The  restrictiveness  of 
the  control  laws  based  on  the  first  type  of  assumptions  is  compensated  for  by  their  simplicity.  The 
price  paid  for  wider  applicability  of  the  control  laws  based  on  the  second  type  of  assumptions  is  an 
overparametrization,  i.e.,  the  need  to  use  the  estimates  of  not  only  the  actual  parameters,  but  also  of 
fictitious  parameters  representing  products  and  powers  of  the  actual  ones. 

This  report  extends  the  results  of  Taylor  et  al.  (1989)  by  introducing  less  n.  *'  live  matching  condi¬ 
tions.  The  new  adaptive  control  law  avoids  any  overparametrization  and  possesses  the  same  robust¬ 
ness  properties  with  respect  to  unmodeled  dynamics. 

We  consider  the  class  of  nonlinear  plants 

x  =fl(x,a)+Fl(x#)z  +Gt(x,a)u  ,  (1.1) 

pi  =  f2(x#)+F2(.xxi)z  +G2(x,a)u  ,  (1.2) 

where  x  e  RH ,  z  e  Rq  are  the  states  of  the  modeled  and  unmodeled  dynamics,  respectively, 
u  e  Rm  is  the  control  (m  £n),  a  e  Rp  is  the  vector  of  unknown  constant  parameters  and  p.  is  a 
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small  positive  scalar  determining  the  presence  of  the  unmodeled  dynamics.  In  the  assumptions 


below,  and  throughout  the  paper,  By  denotes  an  open  ball  centered  at  y  =0. 

Assumption  1  (Equilibrium).  The  columns  of  Fx,  F2,  Gx,  G2,  and  / 1,  f2,  are  real  smooth  vector 
fields  on  RH  such  that 


/ 1(0^) 
fiiW 


e  Im  { 


Gx(0,a) 
G  2(0,0  ) 


for  all  a  e  Ba  . 


(1-3) 


and  hence,  for  every  a  e  Ba  there  exists  a  control  which  places  the  equilibrium  of  the  system 
(1.1H1.2)  at x  =0,z  =0. 


Assumption  2  (Unmodeled  dynamics).  The  unmodeled  dynamics  are  asymptotically  stable  for  all 
fixed  values  ofxefl,  and  a  e  Ba,  that  is,  there  exists  a  constant  o> 0  such  that 

ReX{F2(x/i)}  £  -c<0  .  (1.4) 

■ 

There  may  be  several  reasons  for  not  including  these  dynamics  in  the  design  model  :  our  inability 
to  accurately  model  them;  our  intention  to  simplify  the  model  and  to  reduce  the  controller  complex¬ 
ity;  the  fact  that  their  states  are  not  available  for  feedback.  While  the  smallness  of  p  implies  that 
the  unmodeled  dynamics  are  fast,  a  singular  perturbation  from  p  >  0  to  p = 0,  that  is. 


0  =  f2ixxi)+F2(x,a)z+G2(x,a)u  ,  (1.5) 

results  in  the  reduced-order  model 

m 

x  =  f (xxi)+G(x,a)u  =  f(xji)  +  ’Zgi{x,a)ui  .  (1.6) 

;=i 

where  f  and  G  are  defined  by 

fix, a)  =  f  Xix,a)-F  x{x,a)F2x  ix,a)f  2{x  ,a)  ,  (1.7) 

G(x,a)  =  Gl(x,a)-Fx(x,a)Fi'(x.d)G2(x,a)  .  (1.8) 


The  columns  of  G  (x  & )  are  assumed  to  be  linearly  independent  for  all  x  e  Bx ,  a  e  Ba .  Note  that 
F{1  (xxi)  exists  by  (1.4)  and  that,  by  (1.3), 
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f(0,a)eIm{G(0,a))  for  all  a  e  Ba ..  (1.9) 

Design  tasks.  Given  that  the  x-part  of  the  state  is  available  for  feedback,  our  objective  is  to  find  a 
feedback  control  u  that  will  regulate  the  state  of  the  plant,  i.e.,  that  will  result  in 

x{t),z{t)  bounded,  limx(r)  =  0,  limz(r)  =  0.  (1.10) 

The  three  separate  design  tasks  are: 

(i)  Ideal  design.  When  the  actual  value  of  a  is  known  and  there  are  no  unmodeled  dynamics 
(p=0),  design  a  controller  in  order  to  achieve  regulation  of  (1.6). 

(ii)  Parametric  uncertainty.  Find  a  method  to  cope  with  the  fact  that  the  parameter  vector  a  is 
unknown. 

(iii)  Unmodeled  dynamics.  Establish  that  the  designed  controller  is  robust  in  the  presence  of 
|i>0. 

The  first  task  is  solved  in  Section  2  by  feedback  linearization  of  the  plant  (1.6).  The  second  task  is 
the  topic  of  Section  3,  where  we  develop  an  adaptive  control  law.  The  third  task  is  the  topic  of 
Section  4,  where  the  full-order  closed-loop  system  is  analyzed  and  robustness  bounds  are  computed. 

2.  FEEDBACK  LINEARIZATION 
The  following  assumption  is  crucial  for  our  ability  to  solve  Task  (i). 

Assumption  3  (Feedback  linearization).  For  all  x  e  Bx ,  a  e  Ba  ,  there  exist  a  state 
diffeomorphism  J=<j >(x;a),  with  <|>(0;a)=0,  and  a  feedback  control  u  =a(x',a)+B~l(x-,a)v ,  with 
B(x\a )  a  nonsingular  m  xm  matrix,  that  transform  (1.6)  into  a  linear  system  in  Brunovsky  canoni- 

m 

cal  form  with  controllability  indices  Jk,- ,  j'=l,...,m  (£/t,  ~n) 

/-I 

x  =  Alx+Blv  .  (2.1) 


Necessary  and  sufficient  conditions  for  the  local  validity  of  Assumption  3  are  well  known.  They 
are  also  necessary  (but  not  sufficient)  for  its  global  validity,  which  may  be  difficult  to  verify 
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(Boothby,  1985).  As  shown  in  Jakubczyk  and  Respondek  (1980)  and  Hunt  et  al.  (1983),  Assump¬ 
tion  3  will  hold  locally,  i.e.,  in  a  neighborhood  UQ  of  the  origin  in  Rn,  if  and  only  if  the  distribu¬ 


tions 

G°  =  sp{gl . gM)  ,  (2.2) 

G‘=sp[g , - ,gmMfg\ . adfgm,  .  .  .  ,adjgx . ad}gm]  ,  i=l . n-1  ,  (2.3) 

satisfy  the  following  conditions  for  all  x  e  t/0,  a  e  Ba  : 

GH~i  =  Rn  ,  (2.4) 

C*  is  involutive  and  of  constant  dimension  m,  ,  i= 0 . n~ 2  .  (2.5) 


Then  there  exist  smooth  functions  $,■  of  x,  dependent  on  the  parameter  a,  that  is, 

<b,  (x  \a ),  with  <j>((0;a)=0  and  linearly  independent  differentials  d  <)>,■,  such  that  for  all 
x  e  U0,  a  e  Ba 

ki  =  card  [rjti  ,j  SO  }  ,  ,  r0=m0,  rj  =mJ-mj^  ,  jZ  1 ,  (2.6) 

themxm  matrix  with  <d^itadf  *gj>  as  its  (i,j)-th  element  is  nonsingular,  (2.7) 

<dtyi,adfg[>=b  ,  /=l,...,ro,  i=\,.../n  ,  (2.8) 

(4>i . <t>i . Lf  '<)>,•,  i=l,...,m)  is  a  local  diffeomorphism  which  preserves  the  origin.  (2.9) 

Note,  in  particular,  that  the  controllability  indices  are  required  to  be  the  same  for  all  a  e  Ba ,  and 
that  m0=m  by  the  linear  independence  of  the  columns  of  G.  The  new  coordinates  x  are  then 
defined  as 

x{  =  Lj~l$i(x;a)  ,  j=\,...,ki  ,  j=l . m.  (2.10) 

The  system  (1.6)  is  transformed  into  (2.1)  by  using  the  diffeomorphism  (2.10)  and  the  feedback 
control 

u  a.(x\a)+B~l(x\a)v  ,  (2.11) 

where 

a (x;a)  = -[L/'ih . L^m]T  ,  (2.12) 
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V/vV  •  •  • 

B(x;a)  =  ■  •  (2.13) 

.  . 

The  nonsingularity  of  the  matrix  B(x\a)  for  all  x  e  U0,  a  e  Ba,  is  a  direct  consequence  of  (2.7) 
and  (2.8).  The  resulting  linear  system  of  the  form  (2.1)  consists  of  m  independent  controllable 
subsystems  in  Brunovsky  form,  called  “the  Jj- subsystems.”  For  each  of  these  subsystems  we 
design  a  linear  feedback  control  law 

v,  =  -  X  yfxj/  ,  (2.14) 

so  that  the  &,-th  degree  polynomials  s^+y^s**  '+... +7,^+7/  are  Hurwitz.  This  results  in  exponen¬ 
tially  stable  xt  -subsystems 


(2.15) 


This  design  solves  Task  (i),  that  is,  it  achieves  state  regulation  in  the  case  where  the  parameter  vec¬ 
tor  a  is  known  and  no  unmodeled  dynamics  are  present. 


3.  ADAPTIVE  REGULATION 

For  the  design  under  parametric  uncertainty  (Task  (ii))  we  now  impose  two  restrictions  on  the  class 
of  reduced-order  models  (1.6). 


Assumption  4  (Linear  parametrization).  For  all  x  e  Bx ,  a  e  Ba , 

f(.x,a)=f0(x)+ZfJ(x)aj  , 

8 i (*.<*)  =  g°(x )  +  X  g{(x )aj  , 
im  1 


(3.1) 


i=l,. ...m  , 


(3.2) 


where  /;(x)  and  g{(x )  (j=0,...j>,  i=\,.../n)  are  real  smooth  vector  fields  on  Rn ,  such  that 
/°(0)=0,  f*(! 0)e  G°(0,a),  (note  that  these  are  a  special  case  of  (1.9)),  and 

g\(x,a) . gm{x.a)  are  linearly  independent. 

■ 

Assumption  5  (Extended  matching).  For  all  x  e  Bx  ,  a  e  Ba , 

f'eG1,  ,  (3.3) 

g/eG°,  y=l,...^,t=l . m.  (3.4) 

■ 

This  extended  matching  condition  is  less  restrictive  than  the  following  well-known  condition  (Barm- 
ish  et  al.,  1983,  Marino  and  Nicosia,  1984),  which  was  used  for  adaptive  control  in  Taylor  et  al. 
(1989): 

Assumption  6  (Strict  matching).  For  all  x  e  Bx,  a  e  Ba , 

f  t  e  G°  ,  g{ e  G°  ,  ,  i=l,.../n  .  (3.5) 


The  meaning  of  the  strict  matching  condition  is  that  parametric  uncertainty  is  allowed  to  enter  both 
the  /-vector  and  the  G -matrix  of  (1.6)  only  in  equations  that  also  include  a  control-driven  term. 
Under  the  extended  matching  condition,  the  parametric  uncertainty  in  the  / -vector  of  (1.6)  can  pre¬ 
cede  the  next  control  term  by  one  integrator,  when  the  system  is  transformed  in  the  form  (2.1). 


Example  3.1.  To  illustrate  the  difference  between  the  two  matching  conditions,  let  us  consider  the 
following  simple  but  realistic  dimensionless  second-order  model  of  a  DC-motor: 


=  i  +M(o,a)  , 
dt 


Jj_  di _ 
Tm  dt 


=  —  CO  — X  +ua 


(3.6) 

(3.7) 


where,  in  normalized  units,  co  is  the  motor  speed,  i  is  the  armature  current,  ua  is  the  armature  vol¬ 
tage  (control  input),  and  Tt  and  7m  are  the  electrical  and  mechanical  constants  of  the  motor.  The 
parametric  uncertainty  in  this  case  is  due  to  an  unknown  parameter  a  in  the  torque  characteristic 
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X(co,a)  of  the  motor  load.  The  system  (3.6)-(3.7)  does  not  satisfy  the  strict  matching  condition. 
However,  the  extended  matching  condition  is  satisfied,  since  the  parametric  uncertainty  precedes  the 
control  input  by  only  one  integrator.  The  adaptive  approach  of  Taylor  et  al.  (1989)  cannot  be 
applied  without  first  treating  the  ratio  Tt/Tm  as  a  singular  perturbation  parameter  and  replacing 
(3.6)-(3.7)  by  the  reduced-order  model 


d  to 
dt 


=  -to +ua  +X(a>,a)  , 


(3.8) 


which  clearly  satisfies  Assumption  6.  This  order  reduction  is  not  necessary  for  the  design  approach 
of  the  present  paper. 


Coordinate  transformation.  As  shown  in  Taylor  et  al.  (1989),  under  the  strict  matching  condition 
the  diffeomoiphism  (2.10)  is  independent  of  the  parameters  in  a,  that  is,  <b, (x \a )  can  be  written 
simply  as  <| >,  (x).  However,  under  the  extended  matching  condition  the  diffeomorphism  depends  on 
a  in  the  following  way: 

xj  =  Lfa1 <t>,(;c)  ,  j=l . k  - 1,  i=\,...jn  ,  (3.9) 

x-‘  =  L%  '  <|>,(j<0+  X  ajLfiL%  2<t >,•(*)  ,  i=l,...,m  .  (3.10) 


Since  the  parameter  vector  a  is  unknown,  the  variables  defined  by  (3.9)-(3.10)  are  not  available  for 
feedback.  To  overcome  this  difficulty,  we  replace  the  unknown  a  by  an  estimate  a{t).  Then, 
instead  of  (3.9)-(3.10),  we  use  the  nonlinear  time-varying  change  of  coordinates 


xj  =  L}oX§i(x)  ,  i'=l,... /n  , 

x  *'  =  LfLjo  2<| ),•  (x,d)  =  Ljo  1  $>(x)+  X  dj(t)LfiLkf\ i  20((x)  ,  i=\,...jn  . 


(3.11) 

(3.12) 


If  we  augment  this  change  of  coordinates  by  the  equations  d  =  d,  we  get  a  time-invariant  mapping 
4>  that  maps  ( x,d )  to  ( x,d ).  Since  (3.9)-(3.10)  defines  a  diffeomorphism  for  every  a  e  Ba  and 
(3.12)  depends  only  linearly  on  d,  the  mapping  is  one-to-one,  onto  and  C°°.  The  Inverse  Function 
Theorem  (cf.  Boothby,  1986)  establishes  that  <jj  is  a  diffeomorphism  on  Bx  x  Ba .  Finally,  the  fact 
that  for  every  a  e  Ba  the  diffeomorphism  (3.9)-(3. 10)  preserves  the  origin,  guarantees  that  x=0 
implies  x  =0  and  vice  versa. 


Because  of  the  dependence  of  (3.12)  on  the  parameter  estimates,  the  new  adaptive  control  law  in 
this  paper  is  more  complex  than  the  one  proposed  by  Taylor  et  al.  (1989).  We  now  develop  this 
new  control  law  by  incorporating  not  only  the  parameter  estimate  a,  but  also  its  rate  of  change  a  , 
into  a  “certainty-equivalence”  form  of  the  feedback-linearizing  control  (2.11). 


Certainty  equivalence  control.  Using  Assumption  5,  (2.8),  and  (3.11),  it  can  be  shown  that 


Lfx{  =  LfQxj  ,  ,  i=l,...jn  , 

Lg  x{  =  0,  j=\,...,ki-l  ,  i=\,..jn  ,  r=l . m. 


(3.13) 


Thus,  in  the  new  coordinates  x,  the  plant  (1.6)  consists  of  m  subsystems,  each  described  by  a  set  of 
equations  of  the  form 


x  1 

*  l 


<i>  k.  \  ~  k.  '  mm 

Xi  -  Xi  +w\(x)a  , 

xf‘  =  a{(x)+a2(x)(!+a3(x)a+dra4(x)a+wl,  {x)a 

+  Z  [P  U*  HPi,  U  )d  (x  )a  +a  r  (3^  (X  )a  ]  , 

r-l 


(3.14) 


where  a  =a-d  is  the  parameter  error,  and  the  expressions  for  w\  ,a‘ ,  (3 j7  ,  j= 1 . 4 ,  i ,r=l,...,m 

are  given  in  the  Appendix.  The  design  goal  now  is  to  find  a  control  which  will  render  the  xr 
subsystems  in  (3. 14)  linear  in  the  parameter  error  a  and  make  them  otherwise  independent  of  the 
unknown  parameter  vector  a .  To  achieve  this  goal,  we  calculate  the  control  u  from  the  equation 


B(x,d)u  =  v -w,(x)d  ,  (3.15) 


with  the  obvious  definitions  for  a j  ,y=l,...,4,  W|  and  B{x,d)  (consider  (3.14)  for  j=l,.../n).  Note 
that  B{x&)  can  be  expressed  as 
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This  form  is  directly  obtained  from  (2.13)  when  the  diffeomarphism  (2.10)  is  replaced  by  the  coor¬ 
dinate  transformation  (3.11)-(3.12). 


An  important  feature  of  the  feedback  control  defined  by  (3.15)  is  that  it  depends  not  only  on  the 
parameter  estimate  a,  but  also  on  its  rate  of  change  d  ,  which  appears  in  (3.14)  and  in  (3.15) 
because  of  the  dependence  of  (3.12)  on  a.  The  rate  of  change  d  is  yet  to  be  determined  by  our 
design  of  a  parameter  update  law.  This  update  law  will  assure  that  the  parameter  estimate  a  will 
remain  in  Ba.  This,  in  turn,  will  guarantee  the  invertibility  of  the  matrix  B(x,a)  in  (3.15)  for  the 
reason  mentioned  after  (2.13).  However,  as  d  will  in  general  be  a  function  of  the  control  u,  our 
adaptive  scheme  will  be  implicit,  and,  hence,  the  invertibility  of  B(x,a )  will  no  longer  be  sufficient 
to  guarantee  the  solvability  of  (3.15).  In  the  case  where  a  control  u  satisfying  (3.15)  does  exist,  its 
substitution  in  (3.14)  gives 


x  ■ 

A  l 


-  2  -  *.-l 

-  *.-l  _  k. 

x;‘  =xi'+w\(x)a 

xk<  =  v,-  +w2(x,d,u)a  , 


(3.17) 


with  w2(x,d,u)  as  defined  in  the  Appendix  and  v,-  as  in  (2.14),  but  with  x  replaced  by  x.  A 
further  substitution  of  v,-  into  (3.17)  shows  that  each  of  the  x,- -subsystems  can  be  put  in  the  form  of 
an  error  system 


x,  =  A,x,'  +Wl(x,a,u  )a 


(3.18) 


with  A;  as  defined  in  (2.15).  This  form,  familiar  in  the  linear  adaptive  control  literature  (Narendra 
and  Annaswamy,  1989,  Sastry  and  Bodson,  1989),  will  be  instrumental  in  the  stability  proof.  The 
quantity  multiplying  the  parameter  error  d ,  often  referred  to  as  “the  regressor”,  is 


Wi(x,d,u)  = 


0  .  .  .  0 

o ... .  d 
W  i(x) 

w2(x,d,u ) 


(3.19) 


Parameter  update  law.  The  familiar  form  of  (3.18)  suggests  the  parameter  update  law 
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d  =  -a  =  T  X  Wf(x,d,u)PjXi  ,  (3.20) 

where  T>0  is  the  adaptation- gain  matrix  and  Pi  -Pj> 0  is  chosen  to  satisfy  the  Lyapunov  equa¬ 
tion 

PiAi+AjPi  =  -/*,  ,  i=l,. ...m  .  (3.21) 

Note  that  there  exists  a  control  that  places  the  equilibrium  of  the  system  consisting  of  (1.6)  and 
(3.20)  at  x=0  and  some  d  =  const,  e  Ba.  This  follows  from  (1.9)  and  the  fact  that 
x  =0  o  x  =  0  =5>d  =0. 

Starting  from  (3.20)  and  with  some  algebraic  manipulations,  we  can  show  that  the  parameter  update 
law  is  affine  in  the  control  u ,  that  is, 

d  =rK0(x,d)  +  rKl(x,d)u  ,  (3.22) 

where  (pxl)  and  Kx  (p  xm)  are  defined  in  the  Appendix.  The  matrix  Kx  has  the  property 

Kx{0,d)  =  0  for  all  a  e  Ba  .  (3.23) 

Substituting  (3.22)  into  (3.15)  we  get 

[B(x,d)+Lx(x,d)]u  =v-Lo(;t,d),  (3.24) 

with  the  obvious  definitions  for  L0,  Lx.  From  (3.23)  it  follows  that  Ll(0,d)  =  0  for  all  a  e  Ba,  and, 
since  B(x,d)  is  invertible  for  all  x  e  Bx,  d  e  Ba ,  we  have  the  following  result: 

Proposition  1  (Local  solvability).  For  every  d  e  Ba  there  exists  a  ball  B0  such  that  the  matrix 
B(x,d)+L x{x4)  in  (3.24)  is  invertible  for  all  x  e  B0.  Hence,  in  Ba  there  exists  a  unique  solution 
of  (3.22)  and  (3.24)  (with  v  as  defined  after  (3.17)): 

u  =d(x,d) ,  (3.25) 

<2  =M Kx,d).  (3.26) 

Use  of  this  solution  in  (1.6)  results  in  the  adaptive  closed-loop  system 
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x  =  f(x,a)+G(x,a)u(x,d) , 
d  =y(x4)  . 


The  difficulty  is  that  the  existence  of  such  a  solution  cannot  be  guaranteed  for  all  x  e  Bz ,  a  e  Ba , 
since  in  general  there  exist  values  of  x  and  d  for  which  the  matrix  B(x,a)  +  Lx(x,a)  is  singular. 
One  may  be  able  to  handle  this  singularity  by  a  modification  of  the  update  law.  For  example,  the 
parameter  estimates  can  be  kept  constant  at  their  most  recent  values  while  the  state  is  on  a  manifold 
where  B+L\  is  singular.  For  such  a  modification  to  be  effective,  no  complete  trajectory  of  the 
closed-loop  system  (3.27)  is  to  belong  to  such  a  singularity  manifold,  that  is,  no  singularity  mani¬ 
fold  is  to  contain  an  invariant  set  of  (3.27).  We  stress  that  (3.24)  is  solvable  in  the  whole  region 
bounded  by  the  singularity  manifolds,  and  not  in  some  infinitesimally  small  region,  as  the  term 
“local”  may  suggest.  An  illustration  of  this  fact  is  given  in  Example  3.2  at  the  end  of  this  section. 

Special  cases.  There  are  two  important  special  cases  where  the  global  solvability  of  (3.24)  can  be 
guaranteed.  First,  when  the  strict  matching  condition  (Assumption  6)  is  satisfied.  In  this  case, 
(3.15)  does  not  contain  the  d  -dependent  term,  and  has  a  unique  solution  due  to  the  invertibility  of 
B(x,d).  The  second  special  case  is  when  the  following  assumption  is  satisfied: 

Assumption  7  (Extended  matching  with  parameter-independent  G).  For  all  x  e  Bx ,  a  e  Ba , 

/■'eG1,  p  ,  gj= 0,  ,  i=l . m.  (3.28) 


In  this  case,  the  update  law  (3.20)  is  independent  of  the  control  «,  i.e.,  Ki(;t,d)=0  in  (3.22). 
Hence,  L{(x,d)= 0  and  (3.24)  becomes  uniquely  solvable  due  to  the  invertibility  of  B(x,a ). 

Stability.  The  stability  properties  of  the  closed-loop  system  (3.27)  are  now  established  using  the 
Lyapunov  function 

m 

V(xji)=  'ExlPixi+aTY~la  ,  (3.29) 

i-l 


with  Pi ,  ,  and  T  as  defined  in  (3.21).  The  time  derivative  of  V  along  the  solutions  of 

(3.27)  is 
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V(.x,d)  =  -'£xlxi+2'£xlPiWia+2aTr-la  .  (3.30) 

;-i  i-i 

It  is  now  clear  that  the  update  law  (3.20)  is  chosen  to  guarantee  that 

=  -  S  II II2  =  -  P  INO  ,  (3.31) 

i-i 

and,  hence,  that  x(t) ,  d(t)  are  bounded.  Therefore,  x(t),d(t)  are  bounded,  and  u  is  also 
bounded,  since  the  i=l,...,m,  as  defined  by  (3.25)  and  (2.14),  are  continuous  functions  of  x  and 
d.  By  LaSalle’s  invariance  theorem  (cf.  Hale,  1980)  the  state  of  the  system  (3.27)  will  converge  to 
the  largest  invariant  set  contained  in  the  set  on  which  V(x,a)  =  0,  that  is,  it  =0.  We  conclude  that 
as  t— »<*>,  x(t)— >0.  Since  0,  defined  after  (3.12),  is  a  diffeomoiphism  on  Bx  x  Ba  with 
<J>(0,d )  =  ( 0,d ),  the  fact  that  x(t)— » 0  as  r  — » <»  implies  that  x  (t )  -»  0  as  r  — » 

We  have  thus  proved  the  following  result: 

Theorem  1  (Adaptive  regulation).  Under  Assumptions  1,  3  and  4,  the  state-feedback  control 
(3.25)  and  the  update  law  (3.26)  guarantee  that  the  solution  (x(t),d(t))  of  the  closed-loop  system 
(3.27)  has  the  properties 

lim  x(t)  =  0  , 
x(t),  d(t)-a  bounded, 

for  every  initial  condition  (x0,d0)  that  belongs  to  Q,  where 

(i)  if  Assumption  5  holds,  Cl  is  a  sufficiently  small  neighborhood  of  the  point  (0 ,a)  in  Rn+P , 
or, 

(ii)  if  either  Assumption  6  or  7  holds,  £2  is  the  set 

Q=  {Cc,d):VCc,d)<c}  ,  (3.32) 

with  V(x ,d)  as  defined  in  (3.29)  and  c  the  largest  constant  for  which  the  set  Cl  is  contained  in 


BxxBa. 
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We  examine  now  two  important  special  cases  of  Theorem  1. 

Corollary  1  (Global  adaptive  regulation).  When  Assumptions  1,  3,  4  and  either  6  or  7  hold  for 
all  ( x,a)e  Rn+P ,  the  adaptive  regulation  is  global,  that  is,  Q =RnJfp . 

■ 

Corollary  2  (Parameter  convergence).  When,  in  addition  to  the  Assumptions  of  Theorem  1 ,  the 
following  conditions  are  valid 

(i)  g/=0,  7=1,...,/?,  t=l . m  (parameter-independent  G), 

(ii)  p  <m, 

(iii)  the  vectors  are  linearly  independent, 

the  parameter  estimates  converge  to  their  true  values,  that  is,  the  point  (0,a)  is  an  asymptotically 
stable  equilibrium  of  the  closed-loop  system  (3.27),  with  a  region  of  attraction  Q.  as  defined  in 
Theorem  1. 

Proof.  It  is  sufficient  to  prove  that  the  only  invariant  set  contained  in  the  set  x  =0  is  the  point 

x=0.d  =  0. 

Using  condition  (i),  it  is  readily  shown  that  with  x  =x  =0  the  error  form  (3.18)  reduces  to 

W,(0,a)a=0,  (3.33) 

where 

LjiLfLfl  (j>i(0,<3)  .  .  .  LjpLjLja  <J>i(0,d) 

Ws(0,d)=  ■  (3.34) 

Lf  iLfLjo  "2  <f>m  (0,d )...  LfPLfL{f2  4>m  (0,d ) 

This  is  to  be  compared  with  (3.16)  evaluated  at  x  =0: 
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Lg  JLjLjo  2 <t>i(0,<? )  .  .  .  LgmLfLkfh  2<))1(0,d) 

B(0,d)=  ■  ■  (3.35) 

LgiLfLkff\m(04)  .  •  •  LgmLf-Lkff\m(0,d) 

Condition  (ii),  combined  with  /;(0)e  G°(0)  from  Assumption  4  (with  G°  independent  of  a 
because  of  condition  (i)),  implies  that  the  vectors  /;(0)  can  be  expressed  as  linear  combinations  of 
the  vectors  g,(0) ,  i.e., 

F(0)  =  G(0)D  ,  F  =  [fx  ■  ■  -fp]  .  (3.36) 

with  D  a  constant  m  xp  matrix.  Since  F(0)  is  of  rank  p  by  condition  (iii)  and  G(0 )  is  of  rank 
/n ,  D  is  of  rank  p .  Comparing  (3.34)  with  (3.35)  and  using  (3.36)  we  obtain 

Ws(0,d)  =  B(0,d)D  .  (3.37) 

Since  B(0,d )  is  nonsingular  for  all  d  e  Ba  by  Assumption  3,  and  D  is  of  rank  p ,  (3.33)  and  (3.37) 
imply  that 

a  =0  ,  (3.38) 

and  the  proof  of  the  Corollary  is  complete. 

■ 

Remark  3.1.  Under  the  conditions  of  Corollary  2,  our  adaptive  scheme  achieves  exact  feedback 
linearization  of  the  nonlinear  system  (1.6)  as  t  and  asymptotic  stability  of  the  equilibrium 
x  =  0,  d  =a  of  the  adaptive  system  (3.27). 

Remark  3.2.  In  the  special  case  where  m=p=\  and  Assumption  6  (strict  matching)  is  satisfied, 
condition  (i)  can  be  removed,  that  is,  G  can  be  allowed  to  depend  on  a . 

Example  3.2.  The  system 

ij  =  (10f a)x2,  (3.39) 

•*2  =  ~x2  +(10+a)«  ,  (3.40) 

demonstrates  the  problems  that  can  arise  with  the  solvability  of  (3.24).  In  this  case,  we  have 
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f(x)  = 


(10+a)x2 

_r3 

x2 


,f°(x)  = 


10x2 

r3 

~X2 


. /*(*)  = 


x2 

0 


g(*x)~  [(10fa)]’*0(*)-  [lo]’gl(x)=  [l]’ 


fo  10+a 

f  o  1 

-<10+a)2 

f - 

O 

it 

I _ 

L 10+aJ  ~ 

3(10t-a)x2 

***-  dx 

G°=,p{[i^]}  ,  dim  G°=  1  , 


r  °  i 

-(10+a)2 

LlO+aJ  ’ 

3(10 +a)xi 

}=fl2  for  I  a  I  <  10. 


Thus,  all  the  assumptions  of  Theorem  1  are  satisfied  for  a  e  Ba  =  [a  :  I  a  I  <  10}  and  for  all  x  e  R2. 
Using 


<S>iOO=*i  ,  i./  <t>i==[l  0] 


(10fd)x2 

V3 

~x2 


=  (io +d)x2 , 


the  new  coordinates  are 


Xi=Xi, 

x2  =  (10fd)x2 

When  in  the  transformed  system 

Xi  =  x2+ax2  , 
x2  =  -(l(h-d)x2  +(10f<2)(10+a)«  +d  x2  , 

the  certainty-equivalence  control  is  defined  as 

u  = - - — -[— jcj  —  (10+-<J)j:2— x2  +  (\0+d)x2  ]  ,  (3.43) 

(lOfd)2 

the  system  is  reduced  to  the  error  form  (3.18) 

Xi  =  x2+a  x2  , 

x2  =  -Xi-x2+a(10+d)u  , 

which,  for  a  =0,  is  exponentially  stable.  The  Lyapunov  function 


(3.41) 

(3.42) 
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V(ac^J)  =  °'5]  x 

results  in  the  update  law 

A  =  1.5x1x2  +  (lOh3)[0.5x1  +  (lO4-<!)x2]u  ,  (3.44) 

which,  combined  with  (3.43),  implicitly  defines  the  adaptive  feedback  control  as 

(10fd)[(l(Hd)+x2(0.5A:i  +  (10+d)x2)]u  =  -x,-(l(>*2)x2-1.5xixf  +  0.5(10+-d)x23  .  (3.45) 

For  (3.45)  to  be  solvable,  the  term  multiplying  the  control  u  must  be  nonzero,  which  cannot  be 
achieved  globally.  Indeed,  the  system  (3.39)-(3.40)  satisfies  neither  Assumption  6  nor  Assumption 
7  and  thus  does  not  fall  under  case  (ii),  when  Theorem  1  provides  a  global  result.  However,  Propo¬ 
sition  1  is  satisfied  and  (3.45)  is  solvable  locally.  Suppose  it  is  known  that  a  belongs  to  the  inter¬ 
val  [-8,-2],  and  d0  is  chosen  inside  that  interval.  Then,  if  x, ,  x2  remain  close  to  zero,  the  term 
(10+d)  dominates  the  x  -dependent  term,  thus  guaranteeing  the  solvability  of  (3.45).  Furthermore, 
(3.43)  and  (3.44)  result  in  V(x£)=-  ||x  ||2  ,  thus  yielding  x(r),  d(t)  bounded  and 

lim  x(t)  =  0  ,  lim  x(r)  =  0  . 

i— h»>  »-»«• 

On  the  other  hand,  if  xltx2  are  arbitrary,  (3.45)  may  not  have  a  solution.  In  particular,  if 
d  =-5,  X!  =  15 ,  x2=-2,  the  term  multiplying  u  is  zero.  Nevertheless,  the  region  of  attraction  Q 
defined  in  case  (i)  of  Theorem  I  is  not  some  infinitesimally  small  neighborhood  of  the  origin.  For 
example,  a  subset  of  this  region  is 

{(xxjc2,d):\xx\<2  ,  lx2l  <2  ,  15+4  I  <3}  .  (3.46) 

For  case  (i)  of  Theorem  1 ,  such  (conservative)  estimates  of  the  region  of  attraction  can  always  be 
found. 

4.  ROBUSTNESS  OF  ADAPTIVE  REGULATION 

We  now  return  to  the  full-order  system  (1.1)-(1.2).  The  results  of  Section  2  apply  to  this  system 
when  p.=0.  To  deal  with  the  case  p>0,  we  view  the  z-part  of  the  state  as  consisting  of  two  terms 

z=%(x#,&)+ ri  •  (4.1) 
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Here  ^(x,a,d)  is  the  so-called  “manifold function ” 

^,(x,a,d)  =  -F2l  (x,a)[f2(x,a)  +  G2(x,a)u(x,d)]  ,  (4.2) 

which  is  exactly  equal  to  z  when  jx  =  0,  as  we  can  see  by  comparing  (4.2)  with  (1.5).  The  new 
variable  r\  describes  the  " off-manifold ”  behavior  of  z.  Using  (4.1),  our  plant  (1.1)-(1.2),  along 
with  the  adaptive  controller  defined  by  (3.25)  and  (3.26),  can  be  expressed  as 


x  =  f  (x ,a)+G (x #)u(x ,a)+ F {(x ,a)r\  , 

(4.3) 

pf|  =  F2Cr,a)Ti-p4Ct,<J1<J) , 

(4.4) 

Ci). 

II 

'S? 

a> 

■ — * 

(4.5) 

Our  robustness  result  is: 

Theorem  2  (Robustness).  Suppose  Assumptions  1,  2,  3,  4  and  either  6  or  7  are  satisfied.  Then 
there  exists  a  scalar  p*  >0  such  that  for  all  p  e  [0,p*)  a  Lyapunov  function  V(x,a ,r|)  can  be  found 
which  guarantees  that  (0,0,a )  is  a  stable  equilibrium  of  the  full-order  closed-loop  system  (4.3)-(4.5) 
and  that  the  state  regulation 

lim  *(;)  =  0  ,  limn(r)  =  0,  (4  6) 

is  achieved  for  all  initial  conditions  (x0,T]0,d0)  in  the  set 

C2=  (x,r],d:V(x,d,ri)<c}  ,  (4.7) 

where  c  is  the  largest  constant  for  which  the  set  Q  is  contained  in  Bx  x  Bnx  Ba. 


Proof  (outline).  The  choice  of  the  Lyapunov  function  V  is  made  as  in  Taylor  et  al.  (1989).  The 
requirement  that  the  time  derivative  of  V  for  (4.3)-(4.5)  be  negative  semidefinite  in  Bx  x  B^x  Ba 
is  then  satisfied  by  choice  of  positive  bounding  constants  Cj,  c2  and  c3.  Once  these  constants  are 
determined,  it  is  shown  that  for  Theorem  2  to  hold,  p  can  be  as  large  as 


+ 


1 

CiC2+c3 


(4.8) 


The  exact  form  of  V  and  the  details  of  the  proof  are  given  in  Kanellakopoulos  (1989). 
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Remark  4.1.  In  the  presence  of  unmodeled  dynamics  global  regulation  cannot  be  achieved,  that  is, 
the  results  of  Theorem  2  cannot  be  extended  beyond  Bx  x  Bn  x  Ba  to  Rn+q+p  even  if  its  assump¬ 
tions  hold  globally.  Such  an  extension  would  render  the  constants  Ci.c2.c3  infinitely  large,  thus 
reducing  the  allowed  p-interval  in  (4.8)  to  the  point  0.  The  conflict  is  obvious:  the  larger  we  want 
to  make  the  region  of  attraction,  the  more  we  have  to  restrict  the  allowed  dynamic  uncertainty.  In 
the  limiting  case  p.=0  we  return  to  the  setting  of  Section  2  and  Corollary  1  becomes  valid  again. 

Remark  4.2.  The  nonadaptive-adaptive  comparison  in  Taylor  et  al.  (1989)  applies  mutatis 
mutandis  to  the  adaptive  control  scheme  proposed  here. 


Remark  4.3.  All  the  assumptions  of  Theorem  2,  except  for  Assumptions  1  and  2,  are  required  to 
hold  only  for  the  reduced-order  model  (1.6)  and  not  for  the  full  order  plant  (1.1)-(1.2).  The  benefit 
of  this  is  that  the  assumptions  can  be  satisfied  by  a  broader  class  of  plants. 


Example  4.1.  The  claim  made  in  the  last  remade  is  now  illustrated  by  the  system  consisting  of  the 
DC-motor  of  Example  3.1  driven  by  an  amplifier: 


-^7^-  »  i  +  A,((0,a )  , 
dt 

(4.9) 

Te  di 

—  =  -0)-i  +ua  , 

T„  dt 

(4.10) 

dua 

p—  =  -  ua+k  u  , 

(4.11) 

where  ua  ,k  ,u  ,  are  the  amplifier 
order  system  (4.9)-(4.11)  satisfies 
time  constant  p.  of  the  amplifier  is 
singular  perturbation  from  p.>0  to 


output,  gain  and  input  (control  voltage),  respectively.  The  full- 
none  of  the  Assumptions  5,  6,  7.  However,  observing  that  the 
very  small,  we  can  treat  (<o,0  as  x  and  ua  as  z.  Performing  the 
p.=0  we  obtain  the  reduced-order  model 


=  i  +X(co,a)  . 

at 

di  , 


(4.12) 

(4.13) 


which  clearly  satisfies  Assumption  7.  If  the  load  torque  X(co,a)  is  linear  in  a,  the  system  (4.12)- 


19 


(4.13)  satisfies  all  the  assumptions  of  Theorem  2,  and  can  thus  be  regulated  with  our  approach.  In 
order  to  apply  the  approach  of  Taylor  et  al.  (1989),  one  would  have  to  treat  co  as  x  and  (i  ,ua )  as  z , 
neglecting  the  dynamics  of  both  (4.10)  and  (4.11).  This  would,  however,  result  in  either  more  res¬ 
trictive  robustness  bounds  or  a  smaller  region  of  attraction. 

Let  us  now  return  to  the  exact  setting  of  Example  3.1  by  removing  the  amplifier  (4.1 1).  We  see  that 
the  system  (3.6)-(3.7)  falls  under  case  (ii)  of  Theorem  1,  and  thus  can  be  globally  regulated  using 
the  adaptive  controller  (3.25)-(3.26).  In  contrast,  the  approach  of  Taylor  et  al.  (1989)  provides  only 
a  restricted  region  of  attraction  due  to  the  presence  of  unmodeled  dynamics.  Furthermore,  this 
region  reduces  in  size  as  'he  ratio  Tt  /Tm  increases. 

5.  CONCLUSIONS 

In  the  direct  adaptive  regulation  scheme  proposed  in  this  paper,  the  number  of  parameter  estimates 
is  equal  to  the  number  of  unknown  parameters.  This  is  not  the  case  with  the  general  direct  adaptive 
scheme  of  Sastry  and  Isidori  (1989),  which  requires  that  the  number  of  parameter  estimates  be  aug¬ 
mented.  The  price  paid  for  the  simplification  achieved  here  is  a  restriction  on  the  class  of  plants. 
This  restriction  (Assumption  5)  is  less  severe  than  the  strict  matching  condition  (Assumption  6) 
used  in  Taylor  et  al.  (1989).  It  has  not  yet  been  determined  whether  a  further  relaxation  of  the 
matching  condition  is  possible  without  an  increase  in  the  parameter  estimates. 

Instead  of  broadening  the  class  of  plants  by  a  more  relaxed  matching  condition,  a  similar  effect  is 
achieved  here  by  allowing  the  presence  of  unmodeled  dynamics.  The  same  restriction  is  imposed 
only  on  a  reduced-order  model,  while  the  adaptive  regulation  is  still  guaranteed  for  the  actual  plant, 
provided  that  the  unmodeled  dynamics  are  sufficiently  fast  in  the  sense  of  singular  perturbations. 
As  illustrated  by  the  DC-motor  example,  nonlinear  plants  with  such  unmodeled  dynamics  are  com¬ 
mon  in  practice. 
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APPENDIX 

The  terms  appearing  in  (3.14)  and  (3.17)  are  defined  as  follows: 


a{(x)  =  Lpoto  , 

(A.l) 

P /,(*)  *  LgoL/oX  fo  , 

(A. 2) 

are  scalars. 

it. -2 

w  jj  =  Lp  Lj  o  <J)/  . 

(A. 3) 

4.-2 

®2j(^)  =  lyo  4*i  » 

(A. 4) 

*.-2 

P2 ,r  jOt)  =  LjQ  <(>,  , 

(A. 5) 

*.-l 

0 4j(x)  =  Lfj Lfo  <t»,-  , 

(A. 6) 

p3/j(-*)  =  LgjLfO  4>i  > 

(A.7) 

all  for  J=l,...,p,  are  the  j-th  elements  of  the  1  x  p  row  vectors  wj  (jc),  a^(x), 

Pir.  ajco.  PvCO. 

r=\,...jn,  respectively. 

k.-2 

=  Lft  Lp  L/.‘o  (j),  , 

(A.  8) 

~  Lg‘r  Lfi  L/o  <1 >;  , 

(A.9) 

both  for  l,j=\,...,p,  are  the  (/,y)-th  elements  of  the  p  x p  matrices  a^x), 

P4/U),  r=\,...,m , 

respectively,  and 

w‘2(x,d,u)  =  a3(a:)  +  d7'a4(x)+  'Z[$i,(x)  +  dT$\tr(x)]ur 

r=l 

(A.  10) 

=  Ajr(x,a)  +  uTA?(x,d)  , 

A{(x,d)  =  [LfiXi‘,  .  .  .  ,LfPx-i]T 

(A.  11) 

Ai(x,d)  =  [GlTVxi  *' . G'rVaXiV  , 

(A.  12) 

=  [g{ . g£]  .  y=0,...,p  . 

(A.  13) 

Finally,  the  terms  in  (3.22)  are  defined  as 
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K  0(x,d) 


m  0  ...  0 

L  ,  w? (x) ,  A{(x ,d)  Pi'xi  . 

[  0  .  .  .  0 


Kx(x,d)  =  X(pi*i\  Aj(*4)  . 

<= i 

where  (•)*.  denotes  the  Jfc,  -th  element  of  the  ki  x  1  vector. 


(A.  14) 

(A.  15) 


